A physical interpretation of the cubic map coefficients describing the electron cloud 

evolution 
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The Electron Cloud (ecloud), an undesirable physical phenomena in the accelerators, develops 
quickly as photons striking the vacuum chamberwall knock out electrons that are then accelerated 
by the beam, gain energy, and strike the chamber again, producing more electrons. The interaction 
between the electron cloud and a beam leads to the electron cloud effects such as single- and multi- 
bunch instability, tune shift, increase of pressure and particularly can limit the ability of recently 
build or planned accelerators to reach their design parameters. We report a principal results about 
the analytical study to understanding a such dynamics of electrons. 
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I. INTRODUCTION 

The generation of a quasi-stationary electron cloud in- 
side the beam pipe through beam-induced mult ip acting 
has become an area of intensive study. The analysis per- 
formed so far was based on very heavy computer sim- 
ulations (ECLOUD [![) taking into account photoelec- 
tron production, secondary electron emission, electron 
dynamics, and space charge effects providing a very de- 
tailed description of the electron cloud evolution. 

In has been shown that, for the typical parameters 
of Relativistic Heavy Ion Collider (RHIC), the evolution 
of the electron cloud density can be followed from bunch 
to bunch introducing a " cubic map" of the form: 



*m+l 
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(1) 



where n m is the average electron cloud density after 
TO-th passage of bunch, n m+ i is the one after (m + 1)- 
th passage. The coefficients a, /3, 7 are the parameters 
extrapolated from simulations, and are functions of the 
beam parameters and of the beam pipe features. The 
average longitudinal electron density as function of time 
grows exponentially until the space charge due to the 
electrons themselves produces a saturation level. Once 
the saturation level is reached the average electron den- 
sity does not change significantly. The final decay corre- 
sponds to the succession of the empty bunches. 

A such map approach has been proved, by numeri- 
cal simulations, reliable also for Large Hadron Collider 
(LHC) [3]. The most important outcomes of map for- 
malism for LHC (and generally for any accelerator) are 
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summarizable as follows. In Fig. Q] one can see that the 
bunch-to-bunch evolution contains enough information 
about the build-up or the decay time, although the de- 
tails of the line electron density oscillation between two 
bunches are lost. 
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FIG. 1: Time evolution of the electron density (green line) 
computed with ECLOUD. The black dots mark the average 
electron density between two consecutive bunches. 

Fig. [2] shows the behavior of the average electron den- 
sity n m +i, after the passage of (m+l)-th bunch, as func- 
tion of the average electron density for different 
bunch intensities. 

The points in Fig. [5] show the average electron 
cloud density between two bunches using results from 
ECLOUD, the lines are cubic fits to these points. The 
physical dynamics of electron cloud is explained as fol- 
lows: starting with a small initial linear electron density, 
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FIG. 2: Average longitudinal electron density after the pas- 
sage of bunch as a function of longitudinal electron density 
rim+i before the passage of m-th bunch, for different bunch 
intensities (green N b = 8 * 10 10 , blue N b = 16 * 10 10 ). The 
lines correspond to cubic fits applied to the average bunch 
to bunch points. The red line corresponds to the identity 
map n m +i = n m . Points above this line describe the ini- 
tial growth and saturation of the bunch-to-bunch evolution of 
the line electron density, those below describe the decay. The 
black line represents the cubic fit to the points corresponding 
to the first empty bunches. 



after some bunches the density takes off and reaches the 
corresponding saturation line (n m +i = n m , red line) 
where the space charge effects due to the electrons in the 
cloud itself take place. In this situation, all the points 
(corresponding to the passage of full bunches) are in the 
same spot. The justification of the three terms in eq. 
([I]) is explained as a consequence of the linear growth 
(this term has to be larger than unity in case of electron 
cloud formation), a parabolic decay due to space charge 
effects (this term has to be negative to give concavity to 
the curve rim+i vs n m ), and a cubic term corresponding 
to perturbations. Neglecting the point corresponding to 
the electron cloud density after the first empty bunch, 
the longitudinal electron density follows a similar decay 
independently of the initial value of the saturated line 
electron density. The points corresponding to the first 
empty bunches coming from different saturation values 
lie on a general curve (black curve in Fig. [2]). Thus the 
electron density build up for a given bunch intensity is 
determined by a cubic form, while decay is described by 
two different cubic forms, one corresponding to the first 
empty bunch, and a second to the rest. 

Even though the behavior of the map coefficients from 
simulations has been compared with the experiments al- 
ways this is not well understood and the determination 
of their values is purely empirical. An analytical expres- 



sion for coefficient a (linear coefficient) in the case of 
drift space has been found by understanding the dynam- 
ics that governs weak cloud behavior Q. 

Let us consider iV e /, m quasi-stationary electrons 
gaussian-like distributed in the transverse cross-section 
of the beam pipe. The (m + I)-th bunch accelerates the 
Nei.m electrons initially at rest to an energy £ g . After the 
first wall collision two new jets are created: the backscat- 
tered one with energy £ g and proportional to 5 r (£ g ), and 
the "true secondaries" (with energy £ ~ 5eV) propor- 
tional to S t (£o). The functions S's (or SEY: Secondary 
emission yield) give the ratio of emitted secondary elec- 
trons per incident electron. The sum of all these jets 
becomes the number of surviving electrons Nei >m +i (see 
section Hn]) . and we have 

a = **2±i (2 ) 

In this paper we derive, by assuming the previous the- 
oretical expression for linear coefficient (a) 2], a sim- 
ple approximate formula for the quadratic coefficient (/3), 
which determines the saturation of the cloud due to space 
charge, in the electron cloud density map, under the as- 
sumptions of round chambers and free-field motion of 
the electrons in the cloud. The coefficient depends on 
the bunch parameters, and can be simply deduced from 
ecloud simulation codes modelling the involved physics in 
full detail. Results are compared with simulations for a 
wide range of parameters governing the evolution of the 
electron cloud. 

In the section [TT| we calculate the electronic density of 
saturation by imposing a gaussian-like distribution for 
the space charge and requiring the presence of energy 
barrier near to wall of chamber. In the section IIIII we 
report the calculus of linear coefficient and we compute 
the calculus of quadratic term. Once calculated satura- 
tion we pass to estimate theoretically the coefficient (3. 
We conclude (section IIVI) with the comparison with re- 
spect to outcomes of numerical simulations obtained us- 
ing ECLOUD [T|. In the table U we report all parameters 
used for the the calculations. 



II. STEADY-STATE: ELECTRONIC DENSITY 
OF SATURATION 

In the chamber we have two groups of electrons be- 
longing to cloud: primary photo-electrons generated by 
the synchrotron radiation photons and secondary elec- 
trons generated by the beam induced multi-pactoring. 
Electrons in the first group generated at the beam pipe 
all with the radius b interact with the parent bunch and 
accelerated (by a short bunch) to the velocity v/c — 
2N},r e /b, where r e is the classical electron radius and Nb 
is the efficacy value of bunch population. Since we are 
in presence of trains of finite length beams to consider 
a similar analysis to costing hypothesis it should replace 



3 



Parameter 


Quantity 


Unit 


Value 


Beam pipe radius 


b 


m 


.045 


Beam size 


a 


m 


.002 


Bunch spacing 


Sb 


m 


1.2 


Bunch length 


h 


m 


.013 


Energy for S max 


&Q,max 


eV 


300 


- 


£r 


eV 


60 


Particles per bunch 


N b 


10 10 


4^9 


SEY (max) 


^max 




1.7 


SEY (£ -> 0) 


So 




.7 


SEY (£ -> oo) 


Sac 




.15 




c 




1.83 



TABLE I: Input parameters for analytical estimate and 
ECLOUD simulations. 



the value of bunch population with its spatial or temporal 
average. In fact we have 



Nh 



-N h 



Sb 



(3) 



where Sb is the bunch spacing and h is the length of 
bunch. Electrons in the second group, generally, miss 
the parent bunch and move from the beam pipe wall with 
the velocity v/c = \J 2£q /mc 2 until the next bunch ar- 
rives. The velocity is defined by the average energy £q 
of the secondary electrons and, at high Nb, is smaller 
than velocity of the first group. The process of the cloud 
formation depends, respectively, on two parameters: 



2N b r e h 
b 2 



2£n 



mc 



(4) 



(5) 



These parameters are the distance (in units of b) passed 
by electrons of each group before the next bunch arrives. 
At low currents, k << 1, electron interact with many 
bunches before it reaches the opposite wall. In the op- 
posite extreme case, k > 2, all electrons go wall to 
wall in one bunch spacing. The transition to the sec- 
ond regime can be expected , therefore, k ~ 1 where the 
cloud is quite different than it is at low currents. For 
k > 1 , secondary electrons are confined within the layer 
£ < rjb < 1 at he wall and are wiped out of the re- 
gion < rjb < £ close to the beam by each passing 
bunch. This makes the range of parameters (k > 1 and 
2 — k < £ < 1) quite desirable to suppress the adverse 
effects of the e-cloud on the beam dynamics 

The condition of neutrality implies that secondary elec- 
trons remain in the cloud for a time long enough to af- 
fect the secondary electrons generated by the following 
bunches. In other words, the condition of neutrality and 



the quasi-steady equilibrium distribution of the electron 
cloud are justified only for small k. It is not the case 
at the high currents. In this case, all primary photo- 
electrons disappear just in one pass. The secondary elec- 
trons are produced with low energy £q and are locked up 
at the wall. The density of the secondary electrons grows 
until the space-charge potential of the secondary elec- 
trons is lower than £$. The saturation condition can be 
obtained by requiring that the potential barrier is greater 
than electron energy in the point rjb = 1 — £ 



eV(l-0 ~ So 



(6) 



where V is the electric potential generated by the bunch 
and electron cloud and — e is the electron charge. Here 
it needs to calculate the electric potential by assuming 
some model for the electronic density. 

Our system is composed by a chamber with radius 6, 
a bunch with radius a and length h, an electron cloud 
with density p. Let us consider a electron cylindric dis- 
tribution with a radial gaussian density centered in ro as 
follows 



p(r) 



Poe 



where po is fixed by the condition 



2nh / p{r)rdr — —N e i e 



(7) 



(8) 



where N e i is the total number of electrons in the volume 
irh(b 2 — a 2 ). The electric field in the chamber is 



E = 



N h e 1 N e ,e 
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(9) 



The electric potential V(r) defined by the condition 
V{b) = is 



V(r) 



E-dl 



N b e r _ N e i e X- 
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(10) 
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where F(x) 
g = N b /N el , 
fo = r /b, a 



= /- e 

J a 

VQ 



ydy, G(x) = Si ^dy, 
and x — r/b, a — a/b, 

b 



2neQh 

a/b. We note that in the case a >> 
(or er >> 1) and r — we reobtain the uniform elec- 
tron cloud and with a — > we must neglect the radial 
dimension of bunch with respect to one of electron cloud. 
In fact in this case for eq. (p~0|) we would have 



V ud (r) 



■V 



glrix 



(11) 



Obviously the potentials depend on <?, the ratio of the 
densities of the beam and of the cloud averaged over the 
beam pipe cross-section. In Fig. [3] we report the spatial 
behavior of two potentials. The potential (1111) has min- 
imum at r = r m = b^fg and is monotonic for g > 1 
within the beam pipe. For g < 1 it has minimum at the 
distance r m < b, and the condition g = 1 defines the 
maximum density, this is the well known condition of 
the neutrality. The condition formulated in this form is, 
actually, independent of the form of distribution. Similar 
behavior is found also for the gaussian distribution den- 
sity and is compared with respect to previous one (Fig. 



where p is a free parameter. In the case of uniform dis- 
tribution of electron cloud with a similar mathematical 
passages we find the density of saturation 



N, 



n sa t 
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(14) 



In the Fig. [4] we show the behavior of density of satura- 
tion ([IB"]) and (fT4")l . It is obvious in the case of a gaussian 
distribution of cloud we get a estimate of density satu- 
ration greater than that of uniform distribution. In fact, 
the same number of electrons occupies a smaller volume 
(due to the Gaussian distribution). 
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FIG. 3: Plot of V ~ 1 V(x), go]) and dTTJ), in the case of uni- 
form (solid lines) and gaussian (dashed lines) electronic dis- 
tribution for g — -7- 2, a — .04, fo = 0, a = .3. 

By imposing the condition eq. ([6]) we find the critical 
number (saturation condition) of electrons in the cham- 
ber 
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Nh 



(12) 



while the average density of saturation is found by assum- 
ing that electrons are confined in a cylindrical shell with 
inner radius equal to a and external radius to tq + p a 
where p is a free parameter. So 



io-'X 

FIG. 4: Plot of electronic densities of saturation n aat vs Nb, 
(|13j) and (|14[l ). in the case of uniform (solid line) and gaussian 
(dashed lines) electronic distribution for a — .04, fo = 0, 
a = .3 and p = 2 -j- 3 



III. ANALYTICAL DETERMINATION OF 
COEFFICIENTS 

To compute the linear term we can assume that all 
the Nei,jn electrons, distributed in the transverse cross 
section of the beam pipe, gain an energy E g during the 
passage of the bunch m. After the bunch passage, elec- 
trons are accelerated towards the chamber wall and have 
their first wall collision when two new jets are created: 
one with energy E g and N m S r electrons, corresponding to 
backscattered electrons; the second with low energy and 
N m 8t electrons, corresponding to the true secondaries. 
Before (m + l)-th bunch arrives, these two jets perform 
several wall collisions, which in turn create more jets. 
The contribution of all these jets becomes the number of 
surviving electrons, N m+ \. The reflected electrons travel 
across the beam pipe with energy £ g and perform a num- 
ber of collisions with the chamber wall, between two con- 
secutive bunches, that is: 



N, 



n sa t 



I, sat 



nhb 2 [(rQ + p er) 2 — a 2 



(13) 



(15) 



5 



where t so is the time of bunch spacing, 
tf(£ 9 ) 



iry/2£g/m e 



(16) 



is the average flight time and £ g is the energy of electrons 
accelerated by bunch 
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(17) 



Hence, the total number of reflected, high energy elec- 
trons at the passage of m + 1-th bunch is: 



N. 



el\ref 



N, 



el,mOr 



(18) 



The true secondaries electrons produced after the first 
wall collision gives rise to a low energy jet (So). For this 
jet there is no distinction for the true secondaries and 
reflected, since all are produced with the same energy. 
After the i— th wall collision the number of surviving elec- 
trons is: 



N e i\tot = N e ^ m St S r % 1 S to t k ' 
where S to t = S r (£ ) + S t (£o) and 



hj, — 



I'sb 



itf(£o) 



tf{£o) 



(19) 



(20) 



is the number of collisions after the i— th collision. The 
low energy electrons at the passage of to + 1-th bunch is 



where <5oo, S Q , £ r , (, S max , £ Q , max are the parameters 
of models and their values are shown in Table HI The 
coefficient j3 can be found by imposing the saturation 
condition of map (TTJ): 



n sat = an sat + Pn sat 2 
and the map (Q} becomes 
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1 - a 



n S at 



%i+i 



1 - a 
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(25) 



(26) 



In Fig. ((SJ we show the trends of the coefficient 
(1251) as a function of b max for various values of bunch 
population and viceversa. 
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FIG. 5: Analytical prediction of coefficient ft (|25p for values 
Smax = 1.4 t- 2 and p = 2. 
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Finally the total number of survival electrons at to + 1-th 
bunch passage is obtained taking into account both the 
high and low energy contributions: 



N e l,m+1 — N e l. 



S r S + S t S r l+1 S to t ki 



(22) 



and the linear term @ can be written in the form: 



S r b + S t 5tot V 



Stot 1 - S r S 

where n = ^ £$/ £ g . The expressions of 5's used are 
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FIG. 6: Analytical prediction of coefficient j3 (|25p for values 
N b = (4^9)* 10 10 and p = 2. 



5 r (£) = <5oo + (<5o ~ 5^)6-* IV. CONCLUSIONS 

(24) 

C-i+ £q t 2 ax i iors and the outcomes of simulations (ECLOUD code) 



S t (£) = CS £ i In Figs. ([7]) and (|8J) we report the analytical behav- 
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of coefficients a and ft for values in table |TJ It is im- 
portant to note about the quadratic coefficient we had 
a good agreement with simulations because in our the- 
oretical treatment we have a degree of freedom due to 
the choice of charge distribution in the chamber. By re- 
quiring the condition of saturation it needs to choose an 
average density that can be as realistic as possible. A 
similar discourse was not possible for the linear term be- 
cause we started from the result of Iriso & Peggs (||) and 
then we built the working hypothesis for the determina- 
tion of ft. This specification allows us to evaluate the 
apparent discrepancy between the linear and quadratic 
coefficient with respect to their results from ECLOUD 
code. 

The radial profile of electric density n(r), from 0, is 



n(r) 



p{r) N e i, sa t e 2 



1 — 



2nhb 2 F(l) 



(27) 



and one of energy barrier which opposes the electron com- 
ing from the wall, from (If 01) . is 



g(r) = -eV(r) = eV 



g In x 



G(x) 
F(l) 



(28) 



and we report in Fig. ([9]) their behavior. It notes how the 
peak of the energy barrier corresponds to the maximum 
concentration of electrons around the bunch charge. As 



soon as the density of electrons tends to zero near the 
wall of chamber also energy barrier tends to zero. 
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FIG. 7: Comparison between the analytical prediction (solid 
line) of a (Eq. pg]>) and the simulation (points) by ECLOUD 
code. 

We conclude highlighting the main outcomes of this 
paper. The quadratic map coefficient ft is analytically de- 
rived for the evolution of an electron cloud density. The 
expression is in an acceptable agreement when compared 
with results obtained after ECLOUD simulations with- 
out magnetic field. The analysis is useful to determine 
safe regions in parameter space where an accelerator can 
be operated without creating electron clouds. 
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FIG. 8: Comparison of the quadratic coefficient f3 (Eq. (|25|> ) 
derived using ECLOUD simulations (points) and using the 
analysis of previous sections (dashed lines) with p = 2 -5- 3. 
The solid line is the result by assuming an uniform density. 




FIG. 9: Plot of electronic density. (f27)l. W~ 15 n(r) (dashed 
line) and of energy barrier l)28p . 10 15 £(r) (solid line) for 
a = .3, fo = 0. 



